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SWOVSIS 


Cutting stock problem arises in various industries where stock 
sizes are required to be cut into order sizes* The management is inte- 
rested in finding the cutting patterns that will fill an order at mini- 
mum cost* Present work deals with the problems where order sizes and 
stock sizes are rectangular* 


■In the past, work has been done on this problem for guillotine ' 

: , I 

cutting case, and Gilmore and Gomory have suggested a method which gives 

i 

optimal results* But this method of cutting gives large wastage because | 
it permits only limited flexibility for the combinations of demand rec— 


tangles to be cut from a stock plateyV^is thesis considers the problem 


where cutting equipment permits the cutting of all possible combinations j 
of demand rectangles from a stock plate* Such a problem commonly arises | 
in the cutting of thick plates by gas cutting* An attempt has been made 
to develop an heuristic algorithm to solve this problem* This algorithm 
considers all the cutting patterns permitted in two stage guillotine 
cutting* Along wdth these, some more cutting patterns are considered 
which are generated by the heuristic algorithm* 



(ix) 


This method requires solving a number of knapsack problems# 

An efficient branch and search program has been developed for this. 

It is found that this method gives considerable saving in 
scrap as compared to the two stage guillotine cutting method suggested 
by Gilmore and Gomory. But this saving is at the cost of large computer 
time# 



CHA.PTEB i 


INTRODUCTION 

A problem of primary significance to a variety of industries is 
the reduction of trim losses in cutting of materials for the execution of 
business orders. Typically, this problem is found in paper, glass, and 
steel industries. Management receives orders for various shapes and sizes 
of material. These orders are filled from larger sizes available in stock. 
It is generally found that trimming losses due to odd pieces are unavoi- 
dable. The management is interested in generating cutting patterns that 
will fill an order at minimum cost. This is known as cutting stock problem 
in the field of Operations Research. 

In most of the industries, the order sizes are either one or 
two-dimensional. For example, the order may consist of rolls of given 
widths of a material or rectangular pieces of given widths and lengths etc# 
Correspondingly, there are two types of cutting stock problems. 

One-dimensional cutting stock problem arises in cutting of rolls 

of paper, textile industry, cutting of sections or channels, bars, squares, 

flats, etc. in merchant mills of steel plants. The problem is to minimize 

12 

the scrap losses. Gilmore and Gomory * suggested a solution to this 
problem which gives optimal results. 
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Two-dimensional cutting stock problem arises when there is demaxid 
for two diniensional shapes like rectangles, circles etc# The stock sizes 
are generally rectangular# There are different variations of this problem 
depending on the production process# 

Thin sheets are generally cut by two stage guillotine cutting# 

This problem can be treated as equivalent to two one— dimensional cutting 

stock problems and optimal results can be found by the method suggested 

3 

by Gilmore and Gomory « 

In glass industry, the process is continuous and defective glass 
plates are to be cut into demand sizes# The defective plates are coming 
out of a production shop and they are to be cut in the same sequence in 

which they are coming out# The cutting equipment is capable of guillotine 

''4 . ■ 

cutting only# The sizes of the cuts can be varied# Susan. G# Hahn has 
stiggestecl a solution to such a problem# The optimality of the results can 
not be guaranteed because some subjective decisions have to be made# 

Two-dimensional cutting stock problem is of special significance 
to heavy engineering industry^ where thick plates of scarce material are 
required to be cut into smaller order sizes# Generally gas cutting is 
used for the cutting of thick plates. The concept of ”two stage guillo- 
tinn cutting** can be used for thick plates also, and optimal solution for 
this is available# The draw back of this method of cutting is that it 
permits less flexibility for the combinations of order sizes to be cut 
from a stock rectangle, thus resulting in more scrap losses. So there 
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is a possibility of redaction in wastage if an algoriilun is developed idiich 
permits more flexibility. 

At present, tbis problem is usually solved by hand in industries. 
Experienced men regularly achieve good cutting patterns. But this method 
is time consuming and depends on human judgement. Also, it is difficult 
even for au experienced operator to select the "best" combination without 
actually trying some arrangements. 

In this work, an attempt has been made to develop an heuristic 
algorithm which provides more flexibility for the combinations of order 
sizes to be cut from a stock rectangle and also gives the layout of the 
cutting. 

In chapter 2, literature on one and two-dimensional cutting stodc 
problems and knapsack problem has been reviewed. 

Chapter 3, gives the development of the heuristic algorithra for 
solving two-dimensional catting stock problem. An efficient branch and 
search program for solving the knapsack problem has also been developed, 

Besults for some randomly generated problems are reported in 
chapter 4, A comparison has been made with the results obtained for the 
problems by the method suggested by Gilmore and Gomory for two stage 
guillotine cutting case. 



CHAPTEH 2 

LITEMTUBE BEVIM 

2^1 One-Dinensioiaal Cutting Stock Problem ; 

The simplest kind of cutting stock problem is the one-*dii!iensional 
cutting stock problem^ Consider that k standard lengths ,Lg 

if I . 

of a material are available having costs C^, respectively# 

These will be termed as stock lengths# It is assumed that a large number 
of each of the stock lengths is available# An order has been received for 
lengths Ij^, Ig, * 0 ^ 0 , of the material# These lengths will be termed 
as demand lengths! Let the units of length 1^ demanded be It is 

required to find a cutting arrangement to cut the demand lengths from the 
stock lengths such that the total order cost is minimized. 

5 ' ' ' 

This problem was first formulated by Kurt Eisman as a linear 
programming problem# He suggested that all the possible cutting patterns 
should be enume rated. A cutting pattern shows the number of different 
demand sizes that can be cut from a stock size# It can be written as a 
column 'vector having m components. a . . he the i element of j 

, , , 

column. This represents the ntuaber of demand lengths 1^ in j cutting 
pattern. Let a cutting pattern be called as an activity. The total 
nuiaber of activities will denend on the problem. Now the problem becomes: 
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which activities should be chosen eind how raany times an activity should be 
repeated so that the requirements for demand lengths are satisfied at 
minimum cost? It can be expressed as a linear integer programming problem 
as follows; 


Minimize ^ C. x. 


such that 




a, . X. N. * i s 1, 2, . . . , m 
ij J 1 * 


x.>-0, X. isan integer, j = 1, 2, . . .,n 
3 ' 3 

where n is the total number of activities and x. is the number of 

3 

times j activity is repeated. C. is the cost of the stock length cut by j 

3 

activity. 

Due to unavailability of an efficient method for solving integer 


programming problems, Kurt IDisman siiggested that x.'s can be treated as 

3 

continuous variables so that the problem can be solved as a simple linear 

programming problem. The x- '® general be fractional. 

3 

They can be rounded to nearest smaller integers. It will leave some of the 
orders unfilled. 

13 

Gilmore and Gomory ’ did a detailed study of this problem. If 
the integer restriction on the variables is dropped, the problem can be 
expressed as; 


3»1 


C. X. 
3 3 


Minimize 


6 


n 

such that a., x. - x . »• N. , i » 1,2,# 

13 3 n 4 -x 1 ' 


® I 3 * 1,2, . • *«,n 

a 


where x ,• x . , x are surplus variables. 

n+ 1 n+2 n+m 

It has been proved * that if there is a solution to this problem 
in which all the sxa'plus variables are not zeroes, then there exists a solu- 
tion with the same cost in which all the surplus variables are zeroes. So 
the surplus variables can be droj^ed and the problem can be simplified as: 

n 


Minimize 


3=1 


C. X. 
3 3 


n 

such that a., x. « N. , i = 1,2,^,. .j,m 

^ 13 3 1 ’ 

X- 0 > 3 = l>2,.....jn 

■u 

Simplex method can be used to solve this problem^ The main dis- 
advantage of such a' formulation is that even for a small problem, the 
number of columns (n) is too large for a modem computer, 

13 

In order to overcome this difficulty, Gilmore and Gomory' * 
suggested a column generation technique. It is an efficient way of search-^ 
ing the best colxnnn. The best column is found by solving an auxiliary 
problem which is a knapsack problem. 



Let B be the current basis matrix containing m coltnms 
, Q-,,....,Q 1 . Column Q. is a vector containing m elements. 

d TDl_] 1 

be a nonbasic activity a^^, a 2 ,,»..,,a^ with cost C. The vector 


Let 

Q 


y- ■; , — . — , , , ■ , r 

* NOTE: The proof for this has be« 9 D given by Gilmore and Gomory • This 

propertyholds for this problem because of it's special structure. 
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will gire reduction in cost in current iteration if and only if: 


Cg Q>C 




where Cg is a m component row vector consisting of costs correspond- 
ing to basic variables^ Let Cg B ^ ^ then {2^l) can 

be written as: 

V. a- + v^ + . . . . . V a '> C 
112 2 m m ^ 

So a colnion vector Q is to be chosen sncb that 


^1 ^1 ^ ^ * 


« . * •+v a V C # 

m 


( 2 . 2 ) 


and 


“i * h ^2* 


*Ei “niV 

.... (2.3) 

Sixeh a coltiran vector Q can be found by solving following 
knapsack problem 


Maximize 


m. 

i=l 

m 


V. y. 
1 1 


(2-4) 


such that ^ 1^, E ^|lV 2 ,...., kj (2.5) 


Let 


y. is a positive integer 

'n , yt ,• ■ . ■> yt~\ a solution to this problem. If 


this solution does not satisfy conditions i|2,2) and (2.3) then no other 

columns from the original problem will satisfy these conditions. If v^ 

th 

is interpreted as value of i demand piece, then the objective of this 


* Note: In linear programming, for a minimization problem, a column A. 

■ '.*1 

is a candidate for entering into the basis if and only if 
C^^O, Putting z .ssCg this condition becomes 


J 3 


, which is same as condition (2#l)# 
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knapsack problem is to choose a cutting pattern snch that the total value 
of demand sizes in the cutting pattern is maximized# This approach makes 
the linear programming formulation usable for problems of reasonable sizes* 


2 #2 Two-Dimensional Cutting Stock Problem : 

2#2#1 Two Stage Guillotine Cutting : 

In this problem, a large number of stock rectangles, each 

of width W and length L is available ♦ There is an order for amount 

of a demand rectangle of mdth and length 1^# It is required to find 

a cutting arrangement to cut demand rectangles from the stock rectangles 

by using minimum number of stock rectangles^ It can be formulated as a 

linear programming problem as done in one— dimensional case: 

n 

Minimize x 

j-1 


( 2 . 6 ) 


such that 


a. . X. = N. , i = 1 , 2 , 
ij J 1 




fia 


(2.7) 


X. > 0 . 

J -tb. . . 

where a. . is the number of i demand rectangles in 3 cutting 
13 . 

pattern and x. is the number of times 3 ^^ activity (cutting pattem) 
is to be repeated, n Is the total ntunber of cutting patterns - 


Maximize 


( 2 . 8 ) 


The auxiliary problem for column generation is: 
m 

'■ 1 , '*^'1 " , ■ ■ ■ . .. 

,/ : 1«1 ' , , ■ ■ 

such that fy. , y„, is a feasible cutting pattem 

■ ■ I " 1' ■ III J 

. . ' "t-w-'' — 

( 2 . 9 ) 


where m, y. and v. have the same meaning as discussed earlier* 



It was easy to write (2*9) in mathematical form in case 
of one— dimensional cutting stock problem but the same is not true for 
two-dimensional case* 

3 

Gilmore and Gomory have suggested a solution to the auxi- 
liary problem with the assumption of two stage guillotine cutting# In it, 
the plate is first cut into strips along the length of the plate* The 
width of a strip is equal to the width of atleast one of the demand rec- 
tangles# These strips are in turn used for cutting the demand rectangles# 
The problem of finding a cutting pattern for a stock rectangle such that 
total value is maximized can be treated as equivalent to following two 
problems in one-dimensional case: 

(i) Choosixig a cutting pattern for a strip of know width which 
gives maximum value of the strip# 

(ii) Finding out a combination of strips to be cut from the stock 
plate such that total value is maximized# 

2#2#2 Cutting of Defective Plates: 

Consider now the cutting of defective plates# In a conti- 
nuous production process, some of the plates produced may be defective# 

5 

Such a problem commonly arises in glass industries# Susan G# Hahn has 
suggested a solution to this problem# It has been assumed that the cutt- 
ing equipment allows only guillotine cutting# The delivery dates for the 
demand rectangles are given# 
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It is obvious that pieces from a stock plate should be 
cut according to following policies: 

(i) Preference should be given to the demand sizes having early 
delivery dates and higher demands, 

(ii) Preference should be given to the larger area pieces as compared 
to smaller area pieces, because it is more difficult to fit 
larger pieces. At the same time, care should be taken of the 
fact that cutting of a larger area piece from a more defective 
plate gives more wastage. 

5 

Hahn suggests that values should be assigned to the 
demand rectangles based on areas, due dates etc., and d 3 niainic programming 
approach can be used to solve the problem. 

2.2,3 Generalized Two-Dimensional Cutting Stock Probleo ; 

So far, it has been assumed that cutting equipment permits 
only guillotine cutting and demand sizes are rectangular. In the genera^ 
lized version of the problem, the demand sizes can be in any shapes and 
the cutting equipment permits the cutting of all possible combinations of 
demand sizes# Adamo^ricz and AlkvtXko have suggested a two stage algorithm 
to solve this r>roblem# It has been claimed by these authors that the 
results may not he optimal hut generally they are good# 

The first stage of the algorithm is the clustering stage# 

It produces rectangular enclosures for the shapes or for various cluster- 
ings of shapes according to certain policies# It also gives the numbers 
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of different enclosures to be cut. In the second stage, these authors 
have suggested an heuristic algorithm to find the cutting pattern for a 
stock rectangle such that trim losses for the stock plate are minimized 
without overproducing any demand rectangle-. 

The draw back of this method is that while it attempts to 
minimize the trim losses for individual plates, it does not consider the 
effect of one cutting pattern on another. Secondly^ this method is not 
suitable for problems having high demands for sizes; - 

2.2.4 Present Study ; 

In this work, an attempt has been made to develop an 
heuristic algorithm for the generalized version of the two-dimensional 
cutting stock problem. It has been assumed that demand sizes are rectan- 
gular, This algorithm takes into consideration all the cutting patterns 
permitted in two stage guillotine cutting. These cutting patterns are 
slightly modified as discussed in chapter 3. In addition to these, it 
considers some more cutting patterns generated by the heuristic algorithm; 

2,3 Solution of Knapsack Problem ; 

The efficiency of the proposed method depends on the efficiency 
with which knapsack problem could he solved. 

Knapsack problem has been of interest to research workers 
because it arises in various situations like - capital budgeting, network 
reliability and capital investment etc. It often occurs as a subproblem 
in a problem of optimization. 
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7 

Gilmore and Gomory suggested an efficient dynamic programm- 
ing approach for solving this problem^ This method takes into account 
the periodicity of the solution. The draw back of this approach is that 
it requires a large amount of computer memory# 


8 

Greenberg suggested a branch and search technique when each 
variable y^ can take only two values 0 or 1. This method can also be 
applied to integer solution# There will be only a finite number of solu- 
tions so a method based on successive partitioning can be used. The branch 
and search algorithm proceeds by reT)eatedly partitioning a class of feasi- 
ble solutions into smaller and smaller subclasses along a branch till the 
lower bound is reached or a new feasible solution is found. It then back 
tracks and new branches of the tree are developed. Further branching 
from a node is excluded when the lower bound is reached# Algorithm stops 
when all new branches are excluded# 


It requires solving following linear programming problem at each 


node : 


m 

Maximize v. y. 

" 5 ^ 11 

1=1 

■ ■ m 

such that I. y/ <1 L 

1 . 1 

, ' 1*1 ' 


0 <y. < 1 , i = 1,2, 
y. = Z* 


,m 


’i • ^ € ih] 

gives the indices of the variables which have been assigned 


where 


values at the node. Zf is the value of assigned variable y^ at the node. 
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The variables are first arranged according to decreasing 
magnitude of subscripts be reordered such that 

Vi/lj> V 2 /l 2 j> V Then the solution to this problen 

is: 



y * (L - 'S. 1. 

i * JL 

iCr 


S' z* 1. )/i 


The value of the objective function obtained here gives the 
upper bound on all the solutions represented by this node. Two branches 
can be developed from this node, one by assigning = 0 a*id the other 
by assigning y^ = 1 along with the values assigned at the node. 

The advantage of this method is that the Eemory space required 
is not large. This method has been applied to the original knapsack 
problem (2.4, 2.5) in chapter 3. 



CHAPTEa 3 


TWO-DimiSIOI-IAL GUTTI!^ STOCK PBOBLEM-GENERALIZED CASE 


This chapter considers the solution of two-dimensional cutting 
stock problem, vrhen the cutting equipment nermits the cutting of all 
possible layouts of demand rectangles in a stock rectangle , The method 
developed here is an extension of the two stage guillotine cutting method 
of Gilmore and Gomory. The linear programming formulation of the problem 
as sTiggested by Gilmore and Gomory is : 

n 


liinimize 


such that 




= 1, 2, . . .,m 


Since the number of columns is going to be too large for a 
problem of reasonable size, Gilmore and Gomory suggested a column gene- 
ration technique for choosing the best column. It can be found by solv- 
ing the following auxiliary problem: 

m 

maximize v. y. ... ... (3*l) 

. ^ 11 
1=1 


such that j y . , yo>*»**»y "1 is a feasible cutting 

^ 2 “-i pattern (3.2) 

It will be called as problem P. 




15 


The objective of this auxiliary problem is to choose a cutting pattern 
such that the total value of the demand rectangles in the cutting pattern 
is maximized. The method suggested by Gilmore and Gomory for solving 
this auxiliary problem considers only the cutting patterns permitted in 
two stage guillotine cutting. 

An hexiristic algorithm is proposed for solving the auxiliary 
problem in the generalized cutting case. It will not be possible to 
consider all the cutting patterns by this heuristic algorithm also. 

The optimality of the solution will depend on the total number of cutting 
patterns considered. Larger is the number of cutting patterns considered, 
better the results are expected to he. It is possible that some of the 
cutting patterns permitted in two stage guillotine cutting are not gene- 
rated by the heuristic algorithm and vice versa. So to increase the 
total number of columns considered, the column generation method for two 
stage guillotine cutting and the heuristic algorithm for column genera- 
tion will be combined in the final algorithm. The column generation 
method as suggested by Giliao re and Gomory will he slightly modified, 

3»1 Heuristic Algorithm for Column Generation ; 

The difficulty in solving the problem P is due to ihe feasi- 
bility constraint (3.2). One of the possible ways of solving this 
problem could be by writing (3.2) in mathematical form, hut it is not 
easy to do so. Moreover, the aim here is not just to generate a cutting 
pattern hut also to know it*s layout in the stock rectangle. The 



16 


heuristic algorithm also gives information about it*s layout in the stock 
rectangle# 


Since the total area of ihe pieces in a cutting pattern should 
not- exceed ihe area of the stock plate, the auxiliary problem can be 
written as follows: 


Maximize 


m 

^ V. y. 
W-i, 11 
i«l 


(3,3) 


such that 


and 


m 


i=l 


A. y. A 


11 


(3.4) 






corresponds to a feasible 
cutting pattern 


(3.5) 


.th 


where is the area of i ' demand rectangle and A is the total area 

of the stock plate, is an integer which gives the total numiber of 
th 

i demand rectangles in the cutting pattern. Consider the problem (3.3) 
and (3#4) without the constraint of feasible cutting pattern: 


Maximize 


m 


V. y. 
11 


m 

such that ^ A. y. < A 


i«l 


11 


y^> 0, y^ is an integer , 

Let us call this ^’unconstrained” problem as 


Problem P is the 
1 


familiar knapsack problem. Let f 
be a solution to this problem# This solution may not correspond to a 
feasible cutting pattern# However it gives an upper limit on the total 
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value that can he ohtained* 

The objective of the heuristic algorithm is to find a feasible 
cutting pattern j ^ such that the total value of the objec- 

tive function (3#3) is maximized* The first step of the algorithm is to 

find an upner bound solution I yf, I the auxiliary problem 

^ 1 <5 m*j 

P, A better sblxitioa- th£in. tlii« can not be obtained. It attenpts 
to find out a feasible cutting pattern | al", a*, ...... ,a**^ such that a* 


is closer to yf' for each demand rectangle. The number of different 
demand rectangles aps should not exceed their respective demands. 

The heuristic algorithm tries to fit the demand rectangles 

corresponding to the given cutting pattern |~y*, y*, *1 in the stock 

rectangle according to certain policies discussed later in detail. It 
first arranges the demand rectangles in the decreasing order of preference. 
Let it he ^sumed that the demand rectemgles have been renumbered such that 
they are in decreasing order of preference. The heuristic algorithm con- 
siders the rectangles for fitting in this order. The algorithm has been 

th 

designed such that it attempts to fit a demand rectangle, say i , till 

the total number of it becomes equal to 3^ , In other words, 3^* gives an 

bb 

upper hound on the ntsmber of i demand rectangles that can be considered 
for fitting in this stock plate by the heuristic algorithm. 

The algorithm always attempts to fit the first rectangle if 
possible. If it can not be done, then other demand rectangles are con- 
sidered for fitting, A demand rectangle can not be fitted in following 


two situations ; 
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(a) There is not enough space in the stock plate to fit the demand 
rectangle • 

(b ) The total units of the demand rectangle equals it's upper boxm4 
value 

In situation (b), it is possible that some of the demand rec- 
tangles are not being considered for fitting even though there may be 
enough space for them in the stock plate; There should be some mechanism 
for considering such demand rectangles for fitting, otherwise the empty 
space will go waste. It is taken care of in the second iteration of the 
algorithm. In this iteration, for each rectangle is relaxed, that is 
upper bound for each demand rectangle is increased upto it's total demand 
and the algorithm is continued, Thtis the resulting solution ^ 

will represent a feasible cutting pattern. 


3,1 #1 Some Concepts : 

3,1,1a, Identification of a Stock Plate : 

Subscripted variable E will be used to denote 
the length of edges of a plate. For example denotes the 


th 

length of the J edge. The serial number of the edges of a 
plate are defined as : 

J = 2 6 y 



Fig, 3.1: Numbering of Edges of a Stock Plate 
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According to this definition, odd niiaibered edges are parallel 
to the Y - axis and even ntimbered are parallel to the X - axis# 
The plate at any stage of cutting can be identified by the 
values of the subscripted variable E, This system of identi- 
fying the plates is not applicable to the following types of 
situations: 




■ ” ■ (a) ~ 

Fig. 3*2: Plates Not Identifiable in Present System 

In such a situation, the shaded area shown in figures is 

treated as waste, so that the resulting plates become identi- 
fiable . 

A stock rectangle (shown in Fig. 3.1a) will 
be identified by the following values of the subscripted vari- 
able E: 

Ej^ * Width of the stock rectangle = Eg 

Eg = Length of the stock rectangle* E^ 

The widths and lengths of demand rectangles and the stock 
rectangles are defined such that width ^ length. 
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Following information is needed for fitting a 
demand rectangle in the stock plate and to identify the resul- 
ting stock plate: 

(i) Subscripted variable E to identify the stock plate. 

(ii) The edge J of the stock plate along which the rectangle 
is to be fitted^ The rectangle is always fitted along 
an odd numbered edge. 

(iix) The dimensions of the rectangle to be fitted* Let it 
be B X L (3 and L to be explained later). The edge of 
the rectangle corresponding to the dimension B is to be 

placed along J edge of the stock plate. It is known 

that there is enough space in the stock plate to fit 

th 

this rectangle along it*s J edge. 

For illustration, consider that a rectangle is 
to be fitted along 3rd edge of the stock plate shora in Figure 

( 3.3). 



Fig. 3.3: Stock Plate 
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We have J = 3, One of the situations shown in 
Fig. (3.4 ) will occur depending on the dimensions of the stock 
plate and the rectangle to be fitted. The dotted area shows 
the position of the rectangle in stock plate, and the hatched 
area is treated as waste area. 



B • Eg , L < B « Eg , L - B « Eg , L > E^ 

Pig, 3.4: Pitting a Rectangle in Stock Plate 


Subscripted variable E can be easily redefined 


to identify the resulting stock plate. 


3,1,1c. Orientation of Demand Rectangle 


The rectangle of dimensions B x L represents 
either one of the demand rectangles or a group of similar 
demand rectangles, B and L are calculated depending on the 



orientation of the demand rectangles with respect to the stock 
plate. 

Consider that demand rectangle is to he fitted 

along the J edge of the stock plate. There are two possible 

orientations for the demand rectangle. The orientation is 

chosen such that B is nearest possible to (for J > l). 

d 

Following rules will be followed to choose the orientation. 

It is knoxra that the demand rectangle can be fitted in at least 
one of the txvo orientations. 

(i ) If J =1, then width of the demand rectangle is 
placed along first edge of the stock plate. Let 
Wj. denote the width and 1^ the length of I^^ demand 
rectangle, then B = Wj and L = case the demand 

rectangle can not be fitted in this orientation, the 
other orientation is considered in which B = Ij and 
L = Wj . 

(ii) Consider the case when J>1. In this case, the 

' ' "fcll 

objective is to find out a group of I demand rec- 

' ' ; : 

tangles which will occupy maximum of J edge of the 

stock plate. The group of rectangles will again be 

a bigger rectangle. It‘s dimensions are B x L, B and 

L are calculated from following relationships: 
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L = Ij if B = 



(3*6) 


= Wj. otherwise 


(3.7) 


B here gives the maxiimna of J edge of the stock plate 
that can be occupied. The square brackets ^ here 
denote the highest integer rounded down, 

(iii) In the following situation, the criteria for selecting the 
orientation is different. Consider the stock plate shown 
in Fig, (3,5 ), I demand rectangle is to be placed 
along edge of the stock plate. Two possible orienta- 
tions for the demand rectangle are shown in Fig, ( 3»5). 



If (S, - w^) minimum of dimensions of demand rec- 
tangles, then the hatched area shown in Fig, (3»5b) 
will go waste. In such a case if then it is 

desirable to have the orientation as shown in Fig,(3*5a). 
If i?idth of the demand rectangle is kept along the 3rd. 
edge, the waste area will he less# Hence for such cases 
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we have : 


r • 


B 




X 


and L = 


Chooslns: a Proper Rectangle for fitting: 

The ■rectangles are arranged in decreasing order of pre- 

th 

ference. They are considered for fitting in this order. I rec- 
tangle is not considered for fitting if y* = 0 or the total number 
th 

of 1 demand rectangles fitted is equal to 3^. This rule will be 
followed whenever a demand rectangle is to be chosen for fitting in 
a stock plate. 


3.1.2 Siens of Heuristic Algorithm; 

The objective of the heuristic algorithm is to find a feasible 

•f— ^ 

cutting pattern j a*,a*, . . . . . . ,a* such that total value of the demand 

' ■ I 1 ■ , ■ m I 

rectangles in the cutting pattern is maximized. The values of different 
demand rectangles are given. In th« linear programming formulation, the 
values can be found out by following relationship: 



Now the following steps are taken: 

(0) Find out the upper bound solution F yj> /y^ 

unper bound solution can be found by solving the problem 
This cutting pattern may be infeasible. 
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(i) Arrange the demand rectangles in decreasing order of preference. 
The algorithm arranges: the rectangles according to decreasing 
magnitude of subscripts are reordered such that 

^ ^ n/^m* rectangles are considered 

for fitting in this order. A rectangle is chosen for fitting 
according to the criteria discussed in (3.1. Id). 

(ii ) Choose proper demand rectangle I to be fitted on the J edge 

of the stock plate. In the beginning, take J = 1. 

(iii ) Find out if demand rectangle can be fitted along J edge. 

Consider that a rectangle is to be fitted along 3rd edge of the 

stock plate shoYm in Fig. (3.3 ). The demand rectangle can 
not be fitted with it's width along 3rd edge if: 

(a) Wj> Eg 

(b) lj>E4 + Eg 

Similar conditions can be developed for second orientation also, 
Go to step (ix) if it can be fitted, 

(iv) If it can not be fitted along edge, consider other odd 

numbered edges in the clockwise direction. Go to step (ix) 
if it can be fitted along any of these edges. In case the 
rectangle can not be fitted along any of the odd numbered 
edges, it is concluded that the rectangle, can not be fitted 
in the stock plate at this stage. It does not imply that this 
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rectangle can not be fitted in the stock plate at other stages 
afterwards also. 

(v) Go to step (vi ) if this is the last demand rectangle in the 

ordered sequence (l = m), other wise choose next demand rectan- 
gle for fitting. Go to step (iii) if there is such a demand 
rectangle in the ordered sequence which can he considered for 
fitting, 

(vi) Go to step (viii) if y* , i = 1, m. 

(vii) No more demand rectangles can he cut from the stock plate if 

it is in rectangular form. Go to step (x) in such a situation. 
If the stock plate is not in rectangular form, then there is 
possibility of cutting some more pieces from it by modifying 
the stock plate as shown in Fig. (3,6 ): 



Fig, 3.6: Modifying a Stock Plate 

Take J = 1 and choose proper demand rectangle I. Go to step (iii), 

(viii) The first iteration of the heuristic algorithm is complete. The 
total units of a demand rectangle cut so far does not exceed 
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it's upper bound specified in the cutting pattern found by 
solving the problem . We go for second iteration if there 
is possibility of cutting son® more demand rectangles from the 
remaining stock plate, otherwise go to step (x). The upper 
bounds are relaxed in the second iteration. The algorithm 
redefines the upper bounds as: 
y* = IL , i = 1, 2, ,m 

yf =0 if v^c; 0, i = 1,2,. ,ffi 

Take J = 1 and choose proper demand rectangle from the ordered 
sequence. Go to step (iii), 

(ix) Choose proper orientation for the demand rectangle and calcu- 
late B and L as discussed before. Find the layout of this 
rectangle in the stock plate and redefine the subscripted 
variable E to identify the resulting stock plate. 

None of the demand rectangles can be fitted along first 
edge if Ej is less than the minimum of dimensions of the 
demand rectangles. So the stock plate is modified as shown in 
Fig. ( 3.6 ) and the subscripted variable E is redefined. 
Choose next edge of the stock plate and proper demand rectangle. 
Go to step (iii). 

(x) The algorithm col lects all the waste pieces and finds cutting 
arrangements for them. 
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Finally it finds out the numbers of different demand 

rectangles that have been fitted, af denotes the total number 

"fcll 

of i denand rectangles for each i . 

Flow chart for the algorithm is given in Fig. (3,7 ), 

3.2 Modification of Two Stage Guillotine Cutting Pattern ; 


A cutting pattern in two stage guillotine cutting sometimes 
leaves areas which could be further used for cutting demand rectangles. 
For example, a typical cutting pattern for a strip may be; 



Fig. 3,8 ; Two Stage Guillotine Cutting Pattern 

The hatched area shown in above figure is treated as waste, while 
it my be possible to cut some more demand rectangles from this area. It 
is suggested that the heuristic algorithm developed above can be used to 
find out if any mow.- demand rectangles can be cut. Following procedure 

is followed for this: 

(1) Define subscripted variable E to identify the hatched area. 

(2) Put y* = Total demand for i^*^ rectangle, 1 » 
start tron, step (i) of the heuristic algorittas. 

The output of the heuristic algorith. «11 give the uuriber cf 
.Afferent den«.d rectangles that can he cut fro. the hatched area. 




3.3 Solution of Snapgack Problem: 


The catting pattern P y^, 1 in the sten of the 

<c: m ’ * 

}ieit'.risiic algorithm is generated by solving following knapsack problem: 


.th 


kiaicimi.ze 


stich that 


m 


i==l 


m 


i=l 


V. y. 
1 '^i 


A. y. <1 A 
1 


(3.8) 


(3.9) 


0 , jr^ is an integer for each i (S.lO) 

The variables here have been arranged such 

V This problem ctmld be solved by using dynaniie pro— 
grananing or branch and bound approach. Both these methods require a 
large computer memory. In this thesis, branch and search technique has 
been used to solve this problem. 

In this technique a linear programming problem is solved at 
each node. Consider this problen at node q. 

jffl . ' ’ 

Maximize 

m 

such that ; ’C -A- 




1 1 
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*ere (_Ij| . ( 1,2 represent, the set ot 

indicoe of eseigred voriahles at the node and Z* represent, the valne, 

of the assi,gned variables, A node represents a class of feasible solu- 
tions depending on the assigned variables . The branches emanating from 
a Bole divide this class into subclasses which are mutually exclusive 
but exhaustive# For each branch, we have one or more additionally assi- 
gned variables along with these assigned at the node from which the bran- 
ches are emanating # The solution to the linear programming problem is t 

y^ w Z* , i « — 1 

q-1 

y « (a - 3 1. Z*) / A 

-1 jS: " » ^ 1 

y. “ 0 , j = q + 1, q + 

3 

The value of objective fxinction corresponding to this solution 
gives an upner bound on the class of feasible solutions represented by 
this node, 

3.3,1 Branching Rules : 

Let y =s V (a fractional number) be the solution of 

q 

the linear progranuaing problem at a node from which new branches are to 
be developed. The total number of branches emanating from this node 
will be equal to i ^ | "^ 1» and the value of the additional assigned 
variable y for different branches will be 0, 1, 2, ...••> j V I respec— 
tively. Square brackets [" j here denete the highest integer rounded 


down 
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I'irMt, a brancli will be developed from this node with an addi- 
tional » 0. Linear programming problem is solved for the 

rcanltin^r node. Further branching is continued if the solution at the 
nod! ifi fractional and unner bound exceeds the current lower bound. We 
bftcic tr!>ch if the solution is integer or the upuer bound is less than 
lower borm'’, A node having integer solution will be called as a 
t< T'lninnl node, detailed steps are given below: 

(a) ’iut Cj »» 1 , It represents 1st. node. Find solution to the 


linear programming problem at this node-, taking 

empty set. It will give y = A/A . Find lower bound to 

3. X 


solution, 


It becomes v. 

IS 


a/a1.. 


Put Z* = 


Z* is the 


value of assigned variable y^ for a branch coming out from node 1. 
Put q * 2 wid go to step (b) . 


(b) Find ^ , Z| Z*^^. 

variables y^, 5^2'”"’^q-l 


These are the values of the assigned 
"th, 

at q node. Find solution to the 


linear programning problem without integer restriction, 
give ; y ■ (A — ^ A. Zf) / A . Go to step (c) if y 
integer otherwise go to step (d). 


It will 
is an 


(o) Modify the current lower bound if the value of the objective 
function corresponding to this integer solution is more than 
the current lower bound. Go to step (e), 

(d) Find value of the objective function corresponding to the 
fractional solution. This is upper bound for this node. 
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li’ t." "t. p (0 H upp»r bow.d is less than the oaxrent lewr 
h.n.,l, otherwise further branohiug is done from this node. 

!'( . * « (;, q ■ q+l aujj g.(.gp 

^ to ..top it) if all the branches from (q - node have been 

'hvHoped. Otherwise put + 1 and go to step (b). 

(f) ;;fop if q * p. otherwise put q « q - i and go to step (e). 

This twthod ia illustrated by the help of following example: 

ExnmpU*: Kaxlmi'/e 5x, + + 8x„ 

1 W ti 

such that 4Xj + 3Xg + {?Xg<[^9 

The tree diagram obtained for this problem is shown in Fig. (3.9) 
The dotnila imr described below; 

1» Node 1 ! yj ■ 2|, U.B. ■ 11 1, L.B, * 10, Branch from node 1. 

{*0 to noifcf 2. 

2, Node 2 : » 0 ' ^2 * terminal node* Value of objective 

fsincWon » 10 » current lower bound. Back track. 

li. Node 2 ; - 1 , Jg - 5/8, II.B, - lo| . U.B.> current lower bound, 

m bnanohlng is continued* Oo to node 3. 

4. Nodu 3 ! 7,j> . 1, Z| - 0 , yj . s/2, B-B- - 1“- B.B. - current L.B. 

SO back track. 

5. Node 3 ; Z* - 1, Z| - 1 » Tg - 1* It" is terminal node. Value of 

objective function - 1<^ , vhieJi i® current L.B. 

«o mit l-.B. • ■ ■ , , ' 



35 



Pig* 3*9; Tree Diagram for Branch and Search Example. 
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• Back track, because all the branches from node 2 have been 

explored# 

7# Node 2 : Z* = 2 , = l/3, U#B> = current lower bound# 

Continue branching# 

8# Node 3 : Z* = 2, Z| =0 , y^ = l/g , U.B# = 11 ^ current lower 

bound# Continue branching# 

9# Node 4 : Z* = 2, Z* = 0, = 0 , it is terminal node# Value of 

objective function = 10 = current L#B# , back track. 

lO# : Back track, because all branches from node 2 have been 

explored# 

11* • All branches from node 1 have been explored# Stop because 

q = 2 • 

The final solution is the one corresponding to current lower 
bound# It gives :y| = l, y^ = l, y* = l. 

3 #4 Steps of the Final Algorithm for Solving the Cutting Stock Problem : 


These steps are shoTO in Fig. (3#10 ) # The symbols used in the 
flow chart are: 

= Cost vector for the activities in the basis. 

B ■ ■ 

N = Requirement vector consisting of demands of different 



rectangles , N^, . . . . . ,N^ | . 

Solution vector consisting of m elemenits. It gives the 
number of times different activities in the basis are to he 


repeated 
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C ?= Cost vector consisting of n elements when n is the 
total number of activities# 

Other symbols used are the same as explained before# 

The output of the algorithm is: 

(i) Basis, matrix B# 

(ii) Solution vector 

The columns of the basis matrix B give different cutting 
patterns# There will be total m cuttirg patterns# 

Some of the elements of Xg may be fractional* These numbers 

are rounded to nearest lower integers# Bue to this, the demands for some 

of the rectangles may be unfilled# So some new cutting patterns are 

generated to satisfy the demands of these rectangles. The heuristic 

algorithm for column generation can be used to generate these cutting 

patterns# Following steps are ta^en for this: 

th 

(i) Put Yf = unfilled demnd for i rectangle for each i. This is 
ta^^en as an upuer bound solution# 

(ii) Define value of a demand rectangle equal to it*s area, 

V. = A. , i = l,2,####,ia 

(iii) Start from first step of the heuristic algorithm# Do not go 
for second iteration of the algorithm# Find feasible cutting 
pattern j^a* , 

(iv) Redefine the upper bound solution: 

y^zsjT* — a* , i = l,2,###.,m 
\r.ev«y ^C*"!**; ^ 

stop if yl' e 0 for each i , otherwise go to step (iii)* 
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Thus, this algorithm will give different cutting patterns 
and number of times each of them is to be used to fill the order for 
demand rectangles. The unused areas of the stock rectangles are con- 
sidered as waste. 

Some randomly generated problems have been solved by this 
method and the results are reported in chapter 4, 




Fig. 3.i0: Flow Chart for Final Algorithm 



CHAPTER 4 


RESULTS AMD .DISCUSSION 

4«1 Heuristic Algorithm for Generating Cutting Patterns : 

The objective of this algorithm is to generate a cutting pattern 
having maximum value . The values and demands of different demand rec- 
tangles are given. The nujnber of different demand rectangles in the chosen 
cutting pattern should not exceed their respective demands. 

The first step of the algorithm is to find an upper botmd solu- 
tion [ obtains a feasible cutting pattern 

la*, a*, .a* I from this such that al^' is close to y* for each 

L ^ 2’ m_] X 1 

demand rectangle. 

To test this algorithm, a zero waste cutting pattern for a stock 
plate ( Pig 4.l) was chosen. The value of a demand rectangle is taken 
equal to it's area. It's demand is taken equal to it's number in the 
zero waste cutting pattern. These are shown in table(4.l). The size of 
the stock plate is 14 x 10. 

The maximum value of a cutting pattern here is equal to the 
area of the stock plate. So the zero waste cutting pattern can be taken 


as the upper bound solution 
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TABIiE 4.1 


Zero Waste Cutting Pattern 


Demand 

Rectangle 

- : ■ - 

'^dth 

1 1 

Length Value 

I 

{Demand 
» ' 

!"■ ■ 

-f ■ .. ■ . . 

Upper 

bound solution 

1 Solution by 
{Heuristic 

I Algorithm 


1.0 

3.0 

3.0 

5 / 

5 

5 


1.0 

2.5 

2.5 

3 

3 

3 

^3 

1.0 

1.5 

1.5 

1 

1 

1 


1.0 

1.0 

1.0 

2 

2 

0 

“5 

1.0 

2.0 

2.0 

3 

3 

3 


1.0 

4.0 

4.0 

2 

2 

2 


1.5 

2.0 

3.0 

3 

3 

3 


1.5 

3.0 

4.5 

1 

1 

1 


2.0 

3.0 

6.0 

5 

5 

5 

®10 

2.0 

5.0 

10.0 

1 

1 

1 


2.0 

4.C 

8*0 

1 

1 

1 

^"12 

2.0 

2.5 

5.0 

3 

3 

3 

R 

13 

2.0 

2.0 

4.0 

3 

3 

3 

14 

3.0 

3.5 

10.5 

1 

1 

1 

^15 

0.5 

2.0 

1.0 

1 

1 

1 


The cutting pattern obtained by the hexiristic algorithm is showi 
in column 7 of the table(4.l). Corresponding layout is shown in Fig. (4. 2). 
The total value of this cutting pattern is 138 while the upper bound 
value is 140. The cutting pattern generated here ij 


(fAWuR 

■ GWTBMl library 
Act No. A 30008 
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optimal* The scrap loss of 1.43^ is difficult to explain. Since the 
proposed method is heuristic, the results are quite satisfactory from 
practical consideration. The results may not be so good for all the 
problems, but in general it will give good results. It is interesting 
to note here that time taken for generating this cutting pattern is less 
than 2 secimds cn IBM 7044. 

4,2 Algorithm for Solving Cutting Stock Problem; 

The objective of this algorithm is to find a cutting arrangement 
to cut the demand rectangles from the stock rectangles such that the 
total number of stock rectangles required is minimized, that is, the 
wastage is minimized, * 

This algorithm has been tested on some problems reported below. 
The data for the problems have been generated randomly with uniform dis- 
tributions. A comparison has been made with the results obtained for 
the same problems by two stage guillotine cutting method of Gilmore and 
Gomory , 

The data and the results for the problems are shown in tables 
4.2, 4.3 6uid 4.4. 

The scrap losses by the heuristic algorithm are smaller in 
first two problems and same in the third one as compared to the losses 
by Gilmore and Gomory, It is claimed that the scrap losses by the pro- 
posed method will in no case be more. The reason is that this algorithm 



TABIE 4»2 - Problem 1 


A. Dimensions of Stock Rectangle = 63 x 71 

B, Data for the Problem: 


SI. i 

lo# 

Width 

Length 

Demand 

1 


27 

33 

107 

2 


3 

11 

403 

3 


18 

23 

349 

4 


24 

24 

473 

5 


2 

11 

518 

6 


13 

29 

691 

7 


1 

6 

101 

8 


7 

17 

523 


C. Results for Problem 1; 

Heuristic Gilmore and Gomory 

Wastage 6,14^ 9.78^ 

Computer Time 2 mts. 11 sec. 19 sec. 

IBM 7044 



TABtE 4.3 


Problem 2: 


A, Dimensions of Stock Bectnngle = 100,5 x 208 

B . Data for Problem 2 : 


SI .No. 


Length 

Demand 

1 

13.5 

50.6 

2034 

2 

50.6 

85.7 

8317 

3 

63. G 

86.5 

1602 

4 

6.3 

67.7 

4756 

5 

27.2 

68.4 

4121 

6 

53.0 

60.3 

1923 

7 

47.0 

76.7 

4449 

8 

40.1 

86.2 

9392 

9 

19.7 

28.6 

6313 

10 

58.1 

92.7 

3343 

11 

5.7 

76.0 

3242 

12 

35.9 

94.3 

9364 

13 

23.5 

30.2 

7008 

14 

60.2 

88.2 

7773 

15 

24.6 

93.2 

3565 


C , Hesults for Problem 2; 


Heuristic 

Gi Imore and Gomory 

^*'astage 6.1^ 

8.93^ 

Computer Time 4 nts. 57 sec. 

4 mts* 48 sec. 

IBM 7044^^^^ ^ 
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TABLE 4.4 - Problea 3 

A, Dimensions of Stock Plate = 164.8 x 281 

B. Data for the Prohlem: 


SI. No. 

; Width 

1 

\ Length 

1 

1 

J Demand 

1 

1 

46.2 

79 

2226 

2 

80.3 

55,6 

5455 

3 

20.6 

69.1 

2704 


C. Results for Problem 3: 


Heuristic | Gilmore and Gonory 


Wastage 4.16^^ 4.16^ 


Computer Tiioe 2 sec 
IBM 7044 


1 sec 



takes into account all the cutting patterns permitted in t?/o stage 
guillotine cutting# Fot only this, these cutting patterns are further 
modified before entering into the basis# 

In problems 1 and 2, the scrap losses are about 6^ by the pro 
posed method and by Gilmore and Gomory’s approach# Thus there is 
reasonable improve^ment in the scrap saving# However this saving is 
achieved at the cost of more computer time# The results of problem 2 
are interesting, where the computer time taken by both the algorithms 
are almost sane# It is difficult to explain this behaviour# 

The computer time taken by this algorithm will depend on the 
quality of the colums generated# The number of colurons generated is 
going to be very large in case the individual columns generated are 
of poor quality, and hence more cooputer time ♦ Since the column gene 
ration technique is essentially heuristic, it does not guarantee the 
generation of best columns# 

4 #3 Conclusion : 

From the results, it is obvious that there is going to be con 
siderable saving in scrap losses by using the algorithm developed here 
However, it is at the cost of large computer time# 

This method will be very useful for the practical life prob** 
lens where the material to be cut is costly# There, the saving in 
scra.p may more than makeup for the large computer time# 
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^ ^ ^ need for further research on this method so that the 

computer tine is reduced. ' ■ ■ 

4 #4 ConMents on Memory Eequi remen ts for Branch and Search : ■ 

The branch and search technique for solving the knapsack prob— 
leiss needs very less computer memory. It is because: 

(a) It requires storing information for very less number of nodes 
at a time (maximum HH-I nodes, where m is the number of 
variables in the problem). 

(b) The information stored for each node is also very less* While 
making calculations at a node, the information stored at other 
nodes can be utilized. 

4 #5 Suggestions : 

4. 5.1 Improving the Heuristic Algorithm for Colum Generation : 

The 0th step of the algorithm finds upper bound solu- 
tion by solving following knapsack problem: 

n 

.Maximize v.. y. 

■ "i-' 1 ■ 

■; 1*1 . 

' ^ ■ ■ ' ■ 

, S- A.' y. a ; 

. _ , . . . 

1=1 

y^ is a positive integer for each i . 

A better upper bound solution could be obtained if the maximum number 

of each of the demand rectangles that can be fitted in a stock plate 

th 

is found out. Let be the rmxiirtum number of units of i demand 
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rectangle that can be fitted, then the knapsack problem for finding the 
upper bound solution can be modified as 



m 


Maximize 

^ v^ 
i=l 

m 

5^i 

such that 

>■ A. 

y. c A 

r=i ^ 

1 




y, ■4. minimum (U. N. ) | 
•’x — 1 , 1 V 


is a positive integer 


i = 1,2, 


,m 


■fch. 

where H. is the demand for i ' rectangle. 
1 


The additional constraints can be very easily taken into account 
in the branch and search technique. The time taken for finding upper 
bound solution will reduce, because the total number of nodes generated 
Tifi 11 be less . 


4.5.2 Solution -'flien .Demand lieces Have Irregular Shapes: 


Present work is lisaited to the problems where all the 
demand nieces are rectangular. In practical life problems, the demand 
pieces may have various different shapes. One of the possible ways of 
solving this problem could be by finding rectangular enclosures for the 
demand shapes. It can be done by the method suggested by if.. Idamowicz 
and Albano® . Then, the heuristic algorithm proposed here can be used 
to find the cutting arrangements for the rectangular enclosures. 

4.5.3 Finding a Cutting Pattern Havi ng Hot More Than One Unit 
of Any Demand Rectangle ;. 

An interesting problem could be; A number of demand 
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rectangles is given. A cutting patteni for a stock plate is to be found 
out such that the wastage is minimized. The cutting pattern should not 
have more than one unit of any demand rectangle. Such a problem can be 
easily solved by the heuristic for column generation. The only diffe- 
rence will be in finding the upper bound solution in 0th step. In this 
case, the structure of knapsack problem will be the familiar one where 
each variable can take only two values 0 or 1, 
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